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PERFORMANCE EVALUATION OF AN EMERGENCY CALL 
CENTER: TROPICAL POLYNOMIAL SYSTEMS APPLIED TO 

TIMED PETRI NETS 

XAVIER ALLAMIGEON, VIANNEY BCEUF, AND STEPHANE GAUBERT 


Abstract. We analyze a timed Petri net model of an emergency call center 
which processes calls with different levels of priority. The counter variables of 
the Petri net represent the cumulated number of events as a function of time. 
We show that these variables are determined by a piecewise linear dynamical 
system. We also prove that computing the stationary regimes of the associ¬ 
ated fluid dynamics reduces to solving a polynomial system over a tropical 
(min-plus) semifield of germs. This leads to explicit formulae expressing the 
throughput of the fluid system as a piecewise linear function of the resources, 
revealing the existence of different congestion phases. Numerical experiments 
show that the analysis of the fluid dynamics yields a good approximation of 
the real throughput. 


1. Introduction 

Motivations. Emergency call centers must handle complex and diverse help re¬ 
quests, involving different instruction procedures leading to the engagement of 
emergency means. An important issue is the performance evaluation of these cen¬ 
ters. One needs in particular to estimate the dependence of quantities like through¬ 
puts or waiting times with respect to the allocation of resources, like the operators 
answering calls. 

The present work originates from a case study relative to the current project led 
by Prefecture de Police de Paris (PP), involving the Brigade de sapeurs-pompiers 
de Paris (BSPP), of a new organization to handle emergency calls to Police (number 
17), Firemen (number 18), and untyped emergency calls (number 112), in the Paris 
area. In addition to the studies and experimentation already carried out by PP 
and BSPP experts, we aim at developing formal methods, based on mathematical 
models. One would like to derive analytical formulae or performance bounds allow¬ 
ing one to confirm the results of simulation, to identify exceptional situations not 
easily accessible to simulations, and to obtain a general understanding of poten¬ 
tial bottlenecks. In such applications, complex concurrency phenomena (available 
operators must share their time between different types of requests) are arbitrated 
by priority rules. The systems under study are beyond the known exactly solvable 
classes of Markov models, and it is desirable to develop new analytical results. 
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Contributions. We present an algebraic approach which allows to analyze the 
performance of systems involving priorities and modeled by timed Petri nets. Our 
results apply to the class of Petri nets in which the places can be partitioned in 
two categories: the routing in certain places is subject to priority rules, whereas 
the routing at the other places is free choice. 

Counter variables determine the number of brings of the different transitions as a 
function of time. Our first result shows that, for the earliest firing rule, the counter 
variables are the solutions of a piecewise linear dynamical system (Section]^. Then, 
we introduce a fluid approximation in which the counter variables are real valued, 
instead of integer valued. Our main result shows that in the fluid model, the 
stationary regimes are precisely the solutions of a set of lexicographic piecewise 
linear equations, which constitutes a polynomial system over a tropical (min-plus) 
semiheld of germs (Section |^. The latter is a modification of the ordinary tropical 
semiheld. In essence, our main result shows that computing stationary regimes 
reduces to solving tropical polynomial systems. 

Solving tropical polynomial systems is one of the most basic problems of trop¬ 
ical geometry. The latter provides insights on the nature of solutions, as well as 
algorithmic tools. In particular, the tropical approach allows one to determine the 
different congestion phases of the system. 

We apply this approach to the case study of PP and BSPP. We introduce a sim¬ 
plified model of emergency call center (Section]^. This allows us to concentrate on 
the analysis of an essential feature of the organization: the two level emergency pro¬ 
cedure. Operators at level I initially receive the calls, qualify their urgency, handle 
the non urgent ones, and transfer the urgent cases to specialized level 2 operators 
who complete the instruction. We solve the associated system of tropical polyno¬ 
mial equations and arrive at an explicit computation of the different congestion 
phases, depending on the ratio N 2 /N 1 of the numbers of operators of level 2 and 
I (Section . Our analytical results are obtained only for the approximate fluid 
model. However, they are confirmed by simulations in which the original semantics 
of the Petri nets (with integer firings) is respected (Section]^. 

Related work. Our approach hnds its origin in the maxplus modeling of timed 
discrete event systems, introduced by Cohen, Quadrat and Viot and further devel¬ 
oped by Baccelli and Olsder, see |BCOQ9^ IHOvdWOb] for background. The idea 
of using counter variables already appeared in their work. However, the classical 
results only apply to restricted classes of Petri nets, like event graphs, or event 
graphs with weights as, for instance, in recent work by Cottenceau, Hardouin and 
Boimond |CHB14) . The modeling of more general Petri nets by a combination of 
min-plus linear constraints and classical linear constraints was proposed by Cohen, 
Gaubert and Quadrat |CGQ95| |CGQ98| and Libeaut and Loiseau (see |Lib96) '). 
The question of analyzing the behavior of the dynamical systems arising in this 
way was stated in a compendium of open problems in control theory [Pk99]. A 
key discrepancy with the previously developed min-plus algebraic models lies in the 
semantics of the Petri nets. The model of |CGQ95[|CGQ98| requires the routing to 
be based on open loop preselection policies of tokens at places, and it does not allow 
for priority rules. This is remedied in the present work: we show that priority rules 
can be written in a piecewise linear way, leading to a rational tropical dynamics. 

Our approach is inspired by a work of Farhi, Goursat and Quadrat |TWT| , 
who developed a min-plus model for a road traffic network. The idea of modeling 
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priorities by rational min-plus dynamics first appeared there. By comparison, one 
aspect of novelty of the present approach consists in showing that this idea applies 
to a large class of Petri nets, mixing free choice and priority routing, so that its 
scope is not limited to a special class of road traffic models. Moreover, we provide 
a complete proof that these Petri nets follow the rational tropical dynamics, based 
on a precise analysis of the counter variables along an execution trace. Finally, the 
approach of |FGQ11| was developed in the discrete time case. A novelty of the 
present work consists in the treatment of the continuous time. This requires the 
introduction of a symbolic perturbation technique, working with semifield of germs. 
This technique was used in |GG98) for algorithmic purposes. It has been recently 
applied by All amigeon, F ahrenberg, Gaubert, Katz, and Legay to the analysis of 
timed systems [AFG~*~14] . 

The present piecewise affine dynamical systems bear some general resemblance 
with min-plus models of cellular automata, and in particular with the ultradiscrete 
Toda equation studied by Inoue and Iwao [II12j . We believe these aspects are worth 
being further studied. 

The analysis of timed Petri nets is a major question, which has been exten¬ 
sively studied. We refer to |BD91[ lANOll IGRR041 IJJMS11| for a non-exhaustive 
account on the topic, and to [BV06IILRST091 IBJS09| for examples of tools im¬ 
plementing these techniques. An important effort has been devoted to the com¬ 
parison of timed Petri nets with timed automata in terms of expressivity, see for 
instance |BGH~*~0,^ ISrb08| . The approaches developed in the aforementioned works 
aim at checking whether a given specification is satisfied (for instance, reachability, 
or more generally, a property expressed in a certain temporal logic), or at deter¬ 
mining whether two Petri nets are equivalent in the sense of bisimulation. Hence, 
the emphasis is on issues different from the present ones: we focus on the perfor¬ 
mance analysis of timed Petri nets, by determining the asymptotic throughputs of 
transitions. 

Acknowledgments. We thank Regis Reboul, from PP, in charge of the emergency 
call centers project, and Gommandant Stephane Raclot, from BSPP, for the infor¬ 
mation and insights they provided throughout the present work. We are grateful 
to the anonymous reviewers for their comments which helped to improve the pre¬ 
sentation of this paper. 

2. A SIMPLIFIED Petri net model of an emergency call center 

In this section, we describe a call center answering to emergency calls according 
to a two level instruction procedure. In the new organization planned by PP to¬ 
gether with BSPP [RR15| . the emergency calls to the police (number 17), to the 
firemen (18), and untyped emergency calls (European number 112) will be dealt 
with according to a unified procedure, allowing a strong coordination. Another im¬ 
portant feature of this organization is that it involves a two level treatment. In the 
present paper, we limit our attention to the analysis of the two level procedure. We 
defer to a further work the analysis of the unification of the treatment of calls with 
heterogeneous characteristics. Hence, we discuss a simplified model, for academic 
purposes. 

The first level operators filter the calls and assign them to three categories: 
extremely urgent (potentially life threatening situation), urgent (needing further 
instruction), and non urgent (e.g., call for advice). Non-urgent calls are dealt with 
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Figure 1. Simplified Petri net model of the Parisian 17-18-112 
emergency call center (organization in project). Blue arrows do 
not belong to the Petri net and symbolize the entrance and exit of 
calls in the system. 

entirely by level 1 operators. Extremely urgent and urgent calls are passed to 
level 2 operators. An advantage of this procedure lies in robustness considerations. 
In case of events generating bulk calls, the access to level 2 experts is protected 
by the filtering of level 1. This allows for better guarantees of service for the 
extremely urgent calls. Every call qualified as extremely urgent generates a 3- 
way conversation: the level 1 operator stays in line with the calling person when 
the call is passed to the level 2 operator. Such 3-way conversations allow one to 
avoid any loss of information, and were shown to contribute to the quality of the 
procedure |RR15) . Proper dimensioning of resources is needed to make sure that 
the synchronizations between level 1 and level 2 operators created by these 3-way 
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conversations do not create bottlenecks. We focus on the case where the system is 
saturated, that is, there is an infinite queue of calls that have to be handled. We 
want to evaluate the performance of the system, i.e. the throughput of treatment 
of calls by the operators. 

The call center is modeled by the timed Petri net of Figure We describe here 
the net in informal terms, referring the reader to Section for more information 
on Petri nets and the semantics that we adopt. We use the convention that all 
transitions can be fired instantaneously. Holding times are attached to places. 

Let us give the interpretation in terms of places and transitions. The number 
of operators of level 1 and 2 is equal to Ni and N 2 , respectively. The marking in 
places Pi and p 2 , respectively, represents the number of idle operators of level 1 or 
2 at a given time. In particular, the number of tokens initially available in places 
Pi and P 2 is and A^ 2 - The initial marking of other places is zero. A firing of 
transition qi represents the beginning of a treatment of an incoming emergency call 
by a level 1 operator. The arc from place pi to transition qi indicates that every 
call requires one level 1 operator. The routing from transition qi to transitions 
< 72 , 93,94 represents the qualification of a call as extremely urgent, urgent, or non 
urgent (advice). The proportions of these calls are denoted by Text, 7<'ur, and Tadv, 
respectively, so that Text+Tur+Tadv = 1- The proportions are known from historical 
data. The instruction of the call at level 1 is assumed to take a deterministic time 
Text, Tur, Or Tadv, respectively, depending on the type of call. 

After the treatment of a non urgent or urgent call at level 1, the level 1 operator 
is made immediately available to handle a new call. This is represented by the arcs 
leading to place pi from the transitions located below the places with holding times 
Tur and Tadv Before an idle operator of level 2 is assigned to the treatment of an 
urgent call, which is represented by the firing of transition gg, the call is stocked in 
the place located above gg- In contrast, the sequel of the processing of an extremely 
urgent call (transition gg) requires the availability of a level 2 operator (incoming 
arc P 2 —>■ gg) in order to initiate a 3-way conversation. The level 1 operator is 
released only after a time Ttr corresponding to the duration of this conversation. 
This is represented by the arc gy —>■ pi. The double arrow depicted on the arc 
P2 gg means that level 2 operators are assigned to the treatment of extremely 
urgent calls (if any) in priority. The holding times and represent the time 
needed by a level 2 operator to complete the instruction of extremely urgent and 
urgent calls respectively. 

3. Piecewise linear dynamics of timed Petri nets with free choice 

AND priority ROUTING 

3.1. Timed Petri nets: notation and semantics. A timed Petri net consists of 
a set V of places and a set Q of transitions, in which each place p GV is equipped 
with a holding time Tp € M>o as well as an initial marking Mp € N. Given a 
place p G V, we respectively denote by p“ and the sets of input and output 
transitions. Similarly, for all q G Q, the sets of upstream and downstream places 
are denoted by g‘" and g°“‘ respectively. 

The semantics of the timed Petri net which we use in this paper is based on 
the fact that every token entering a place p G V must stay at least Tp time units 
in place p before becoming available for a firing of a downstream transition. More 
formally, a state of the semantics of the Petri net specifies, for each place p G V, 
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the set of tokens located at place p, together with the age of these tokens since they 
have entered place p. In a given state cr, the Petri net can evolve into a new state 
fj' in two different ways: 

(i) either a transition g G Q is fired, which we denote a —^ a'. This occurs 
when every upstream place p contains a token whose age is greater than or equal 
to Tp. The transition is supposed to be instantaneous. A token enters in each 
downstream place, and its age is set to 0; 

(ii) or all the tokens remain at their original places, and their ages are incre¬ 
mented by the same amount of time d S This is denoted a —^ a'. 

In the initial state cr®, all the tokens of the initial marking are supposed to have 
an “infinite” age, so that they are available for firings of downstream transitions 
from the beginning of the execution of the Petri net. The set of relations of the 
form —^ and —^ constitutes a timed transition system which, together with the 
initial state cr°, fully describe the semantics of the Petri net. Note that in this 
semantics, transitions can be fired simultaneously. In particular, a given transition 
can be fired several times at the same moment. Recall that every holding time Tp 
is positive, so that we cannot have any Zeno behavior. 

In this setting, we can write any execution trace of the Petri net as a sequence 
of transitions of the form: 

(1) cr° jLf. 

where ^ 0 and > 0. In other words, we consider traces in which 

we remove all the time-elapsing transitions of duration 0, except the first one, 
and in which time-elapsing transitions are separated by groups of firing transitions 
occurring simultaneously. We say that a transition q is fired at the instant t if there 
is a transition —^ in the trace such that the sum of the durations of the transitions 
of the form —> which occur before in the trace is equal to t. The state of the Petri 
net at the instant t refers to the state of the Petri net appearing in the trace Q 
after all transitions have been hred at the instant t. 

In the rest of the paper, we stick to a stronger variant of the semantics, referred 
to as earliest behavior semantics, in which every transition q is fired at the earliest 
moment possible. More formally, this means that in any state a arising during the 
execution, a place p is allowed to contain a token of age (strictly) greater than Tp 
only if no downstream transition can be fired (i.e. no transition —^ with q € p°'^* can 
be applied to a). The motivation to study the earliest behavior semantics originates 
from our interest for emergency call centers, in which all calls are supposed to be 
handled as soon as possible. 

3.2. Timed Petri nets with free choice and priority routing. In this paper, 
we consider timed Petri nets in which places are free choice, or subject to priorities. 
This class of nets includes our model of emergency call center. Recall that a place 
p G V is said to be free choice if either = 1, or all the downstream transitions 
q G p°'^^ satisfy = {p}. The main property of such a place is the following: if 
one of the downstream transitions is activated (i.e. it can be potentially fired), then 
the other downstream transitions are also activated. A place is subject to priority 
if the available tokens in this place are routed to downstream transitions according 
to a certain priority rule. We denote by Apriority the set of such places. We assume 
that no transition has more than one upstream place subject to priority, that is, for 
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Figure 2. Conflict, synchronization and priority configurations. 


any transition g, the set g™ H T^priority has at most one element. This allows to avoid 
inconsistency between priority rules (e.g. two priority places acting on the same 
transitions in a contradictory way). For the sake of simplicity, we also assume in 
the following that every p € T^priority has precisely two downstream transitions, which 
we respectively denote by p™* and p™*. Then, if both transitions are activated, 
the tokens available in place p are assigned to p™‘ as a priority. Equivalently, 
in the execution trace of the Petri net, we have a —a' only if the transition 

out 

—j-Pt cannot be applied to the state a. We remark that it is possible to handle 
multiple priority levels, up to making the presentation of the subsequent results 
more complicated. 

To summarize, there are three possible place/transition patterns which can occur 
in the timed Petri nets that we consider, see Figure The first two ones involve 
only free choice places, and are referred to as conflict and synchronization patterns 
respectively. We denote by T^confiict the set of free choice places that have at least 
two output transitions, and by Qsync the set of transitions such that every upstream 
place p satisfies |p°“*| = 1. By definition, we have ^’conflict F (Qsync)'" = 0- The third 
configuration in Figure [^depicts a place p subject to priority. In order to distinguish 
p™‘ and p™*, we depict the arc leading to the transition p™* by a double arrow. 
By assumption, the places r p located upstream p™* and p™* are non-priority, 
so that they are free-choice and have only one output transition, as depicted in 
Figure l^c). 

3.3. Piecewise linear representation by counter variables. Since we are in¬ 
terested in estimating the throughput of transitions in a Petri net, we associate 
with any transition g S Q a counter variable Zq from K to N such that Zq{t) rep¬ 
resents the number of firings of transition g that occurred up to time t included. 
Similarly, given a place p G V, we denote Xp{t) the number of tokens that have 
entered place p up to time t included. Note that the tokens initially present in 
place p are counted. More formally, Xp{t) is given by the sum of the initial marking 
Mp and of the numbers of firings of transitions g G p'" which occurred before the 
instant t (included). We extend the counter variables Xp and Zq to K<o by setting: 

(2) Xp{t) = Mp , Zq{t) = 0 , for all t < 0 . 

By construction, the functions Xp and Zq are non-decreasing. Besides, since they 
count tokens up to time t included, they are cddldg functions, which means that 
they are right continuous and have left limits at any point. Given a cadlag function 
/, we denote by f{t~) the left limit at the point t. 
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The goal of this section is to describe the dynamics of timed Petri nets with free 
choice and priority routing by means of a set of piecewise linear equality constraints 
over the counter variables. We provide an informal presentation of these constraints. 
First observe that we necessarily have: 

(3) VpGT’, Xpjt) = Mp + Zqjt ), 

as the initial marking Mp is counted in Xp{t), and any token entering place p before 
the instant t must have been fired from an upstream transition q G p'" before. In 
a similar way, if p G T^confiictj the total number of times the downstream transitions 
have been fired before the instant t is necessarily equal to the number of tokens 
which entered place p before time t — Tp (included). This is due to the fact that if 
a token enters p at the instant s, then it is consumed exactly at the instant s + Tp 
(by definition of the earliest behavior semantics). This yields the identity: 

(4) Vp G "Pconflict , ^ Zq{t) = Xp{t - Tp) . 

q^poMt 

Now consider a transition q G Qsync- The number of times this transition is fired 
at the instant t is given by Zq{t) — Zq{t~). In each upstream place p G q™, the 
number of tokens which are available for firing q is equal to Xp(t — Tp) — Zq(t~). 
Indeed, since place p does not have any other output transition, the total number 
of tokens which have left place q until the instant t equals Zq{t~). By definition of 
the earliest behavior semantics, the number of firings of q at the instant t must be 
exactly equal to the minimum number of tokens available in places p G q™. If we 
denote min(a;, p) by a; A y, we consequently get: 

(5) Wq G Qsync ; ^qi^) ~ l\ ^pi^ '^p) ■ 

peq'’' 

Finally, let us take a place p G T^priority Since the transition p™* has priority over 
the quantity Zp<^t (t'j Zp<^t (t ) must be equal to the minimal number of tokens 
available in the upstream places, including p. For every place r G (p™*)'" distinct 
from p, the number of available tokens is given by Xr{t — Tr) — Zp<^t {t~) (recall that 
p™‘ is the only downstream transition of r). In contrast, the number of tokens 
available for firing in place p is equal to Xp{t — Tp) — {Zp<^t[t~) + Zpout(t~)). We 
deduce that we have: 

(6) Vp G 7^priority r (7) — {xpit Tp) Zpi^tlt )) A Xj-it Tr) . 

r&ipT)"' 

r^p 

The number of tokens from place p which are available for the transition after 
the firings ofp™* is given by Xp{t—Tp) — {zp-^t{t~)+Zpi>ut{t~)) — {zp'^t(^t) — Zp'^t{t~)). 
Hence, we obtain: 

(7) Vp G T^priority r Zpf^t (t) — {xp{t Tp) Zp^ut(t)^ A Xr{t Tr) . 

r^p 

We summarize the previous discussion by the following result: 
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Theorem 1. Given any execution trace of a timed Petri net with free choice and 
priority routing, the counter variables Xp (p G V) and Zq (q G Q) satisfy the 
constraints @-0 for all t ^ 0, together with the initial conditions 0. 

We refer to Appendix [X| for a detailed proof of this statement. Notice that, if 
we do not restrict to the earliest behavior semantics, the constraints 0-0 are 
relaxed to inequalities. 

So far, we have described the dynamics of timed Petri nets in the continuous time 
setting. However, since the Petri net of our case study is a model of a real system 
which is implemented in silico, we need to investigate the dynamics in discrete time 
as well. In more details, assuming that all the quantities Tp are multiple of an 
elementary time step i5 > 0, the discrete-time version of the semantics of the Petri 
net restricts the transitions —^ to the case where d is a multiple of 6. In this case, 
on top of being cadlag, the functions Xp and Zq are constant on any interval of the 
form [kS, {k + l)6) for all k gN. Then, we can verify that the following result holds: 

Proposition 2. In the discrete time semantics, the counter variables Xp and Zq 
satisfy the constraints 0-0 for all t ^ 0, independently of the choice of the 
elementary time step 6. 


In other words, the dynamics in continuous-time is a valid representation of the 
dynamics in discrete time which allows to abstract from the discretization time 
step. We also note that we can refine the constraint given in 0 by replacing the 
left limit Zpout(t~) by an explicit value: 


( 8 ) 


Vp G ^priority , 


' {Xp(t - Tp) - Zpout (t - (5)) 


(t) = < 


A /y Xr{t 

, r^p 

{Xp{t - Tp) - Zpa^t{f)) 


Tr) 


A /y Xr{t — Tr) 

rGlp""*)'" , r^p 


iff e (5N, 


otherwise. 


(Here and below, we denote by dN the set {0, 6,2d,... }.) The system formed by 
the constraints 0-0, 0, 0 is referred to as the 6-discretization of the Petri net 
dynamics. 

The only source of non-determinism in the model that we consider is the routing 
policy in the conflict pattern (Figure 0a)). In the sequel, we assume that the 
tokens are assigned according to a stationary probability distribution. Given a 
free choice place p G 'Pconfiict, we denote by iTqp the probability that an available 
token is assigned to the transition q G In the following, we consider a fluid 

approximation of the dynamics of the system, in which the Xp and Zq are non¬ 
decreasing cadlag functions from K to itself, and the routing policy degenerates in 
sharing the tokens in fractions Equivalently, the fluid dynamics is defined by 
the constraints 0-0 and the following additional constraints: 

(9) Vp G ^conflict , Vq G , Zq{f) = TTqpXp{t - Tp) . 


Note that the latter equation is still valid in the context of discrete time. By 
extension, the system formed by the constraints 0-0, 0-0 is referred to as 
the 5-discretization of the fluid dynamics. 
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3.4. Application to our Petri net model of emergency call center. We 

illustrate Theorem on the Petri net of Figure We point out that in Figure 
we have omitted to specify the holding time of some places. By default, this holding 
time is set to a certain > 0, and is meant to be negligible w.r.t. the other holding 
times. 

For simplicity, we omit the counter variables of the places distinct from pi and 
P 2 - Indeed, each of theses places p has a unique input transition g, and its initial 
marking is 0. Therefore, by definition, we have Xp{t) = Zq(t) for all t, which means 
that Xp can be trivially substituted in the constraints. Similarly, we omit the 
transitions which lead to places pi and p 2 , as their counter variables correspond the 
counter variables of some transitions located upstream and shifted by the holding 
time of the place in between. Finally, we denote by Zi the counter variables of 
transitions qi, and by Xi the counter variables of places Pi. We can verify that the 
fluid dynamics is then given by the following constraints: 

Zi{t) =Xi{t- Te) 

Z2{t) = 7rext2:i(t - Te) 

Zd.{t) = - Te) 

Zi{t) = 7radv-2:i(t - Te) 

Z^{t) = {X2{t - Te) - Ze{t~)) A Z2it - Text) 

Ze{t) = {X2{t - Te) - Z^{t)) A Zzit - Tur “ Te) 

Z7{t) = Z^{t- Ttr) 

Xi{t) = Ai + Z7{t) + Z'iit - Tur) + Ziit - Tadv) 

X2{t) = N 2 + Z7{t - + Zeit - r'J 

They can be simplified into the following system: 

Zi{t) = Ni+ Z5{t - Ttr) + TurZlit - T^r - ^Te) + 7radv2:i(t “ Tadv “ 

Zb{t) = {N 2 + Z^{t - Ttr - Text “ ^e) + Zeit - T^^ - Te) - Ze(t~)) 

(10) A TTextZl (t - Text “ Te) 

Zeit) = {N 2 + Zeit - Ttr - Text “ Te) + Zeit - t(„. - Te) - Zeit)) 

A Tr^rZiit - Tar “ Te) 


which involve the counter variables zi, ze and ze only. These variables correspond 
to the key characteristics of the system. They respectively represent the number of 
calls handled at level 1, and the number of extremely urgent and urgent calls han¬ 
dled at level 2, up to time t. All the other counter variables can be straightforwardly 
obtained from zi, Z5 and ze- 

For the sake of readability, we slightly modify the original holding times Text, 
Tur, ... to incorporate the effect of t^. In more details, we substitute Text, Tur, Tadv, 

TLt and Tur by Text “ T^, Tur “ St^, Tadv “ St^, t'^^ - Tg and T^^r “ Te respectively. 

Then, System (10) simply reads as: 

Zi(t) = Ai -I- Z5(t - Ttr) + 7 rurZl(t - Tur) + T'r.dvZlit - Tadv) 

(11) Zeit) = (A 2 -b Zeit - Ttr - Text) + Zeit - Tur) “ Zeit~ )) A TTextZlit - Text) 

Zeit) = (A ^2 + Zeit - Ttr - Tg,rt) + ^f>it “ T'r) “ Zeit)) A 7rur2:i(t “ Tur) 


This is the system which we consider in the rest of the paper. 
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4. Computing stationary regimes 

We investigate the stationary regimes of the fluid dynamics associated with Petri 
nets with free choice and priority routing. More specifically, our goal is to charac¬ 
terize the non-decreasing cadlag solutions Xp and Zq of the dynamics which behave 
ultimately as affine functions t u + pt (u € M. and p G K^o)- By ultimately, 
we mean that the property holds for t large enough. In this case, the scalar p 
corresponds to the asymptotic throughput of the associated place or transition. 
However, if the functions Xp and Zq are continuous, and a fortiori if they are affine, 
their values at points t and t~ coincide, and then, the effect of the priority rule 
on the dynamics vanishes (see Equation Q). Hence, looking for ultimately affine 
solutions of the continuous time equations might look as an ill-posed problem, if one 
interprets it in a naive way. In contrast, looking for the ultimately affine solutions 
of the (5-discretization of the fluid dynamics is a perfectly well-posed problem. In 
other words, we aim at determining the solutions Xp and Zq of the discrete dynam¬ 
ics which coincide with affine functions at points k5 for all sufficiently large k gN. 
These solutions are referred to as the stationary solutions of the dynamics. As we 
shall prove in Theorem the characterization of these solutions does not depend 
on the value of i5, leading to a proper definition of ultimately affine solutions of the 
continuous time dynamics. 

In order to determine the stationary regimes, we use the notion of germs of affine 
functions. We introduce an equivalence relation ^ over functions from K to itself, 
defined by / ~ 5 if /(t) and g{t) are equal for all t G 5N sufficiently large. A germ 
of function (at point infinity) is an equivalence class of functions with respect to 
the relation ~. For brevity, we refer to the germs of affine functions as affine germs, 
and we denote by (p, u) the germ of the function t i-G u + pt. In this setting, our 
goal is to determine the affine germs of the counter variables of the Petri net in the 
stationary regimes. 

Given two functions / and g of affine germs {p,u) and {p',u') respectively, it is 
easy to show that f{t) ^ g{t) for all sufficiently large t G 5N if, and only if, the 
couple (p, u) is smaller than or equal to {p', u') in the lexicographic order. Moreover, 
the affine germ of the function f + g is simply given by the germ (p -|- p^ u -I- u^), 
which we denote by (p, u) -I- {p' ,u') by abuse of notation. As a consequence, affine 
germs provide an ordered group. Let us add to this group a greatest element T, 
with the convention that T -|- (p, m) = (p, u) -I- T = T. Then, we obtain the tropical 
(min-plus) semiring of affine germs (G, A, -I-), where G is defined as {T} U and 
for all x,y G G, X A y stands for the minimum of x and y in lexicographic order 
(extended to T). Since in G, the addition plays the role of the multiplicative law, 
the additive inversion defined by —(p, u) := (—p, —u) corresponds to a division over 
G. This makes G a semifield, i.e., in loose terms, a structure similar to a field, except 
that the additive law has no inverse. Finally, we can define the multiplication by 
a scalar A € K by A(p, u) := {Xp, Xu). When A S N, this can be understood as an 
exponentiation operation in G. 

Instantiating the functions Xp and Zq by affine asymptotics t Up tpp and 
1 1 —>■ Uq+tpq in the (5-discretization of the fluid dynamics leads to the following coun¬ 
terparts of the constraints (§, (§, 0 and 0 , the variables being now elements 
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of the semifield G of germs: 


(12a) 

ypGV, 

(Pp, Up) 

— (0, Mp) -1- 'y ( (Pq, Uq) 





(12b) 

Vp G (Pconflict ,V(7 G p°'^^ , 

{Pq,Uq) 

~ 'TqpiPp^ '^p PpTp) 

(12c) 

yq G Qsync ) 

{Pq,Uq) 

— /\ {Ppt '^p ~ PpTp) 




(12d) Vp G 1^ priority i : Upc^t'j {Pp Pp^^ : ^p Pp'^p 

A {Pr: '^r Pr'^r) 

, r^p 

Given p G T^priority, the transposition of (|^ (or equivalently ([^) to germs is more 
elaborate due to the occurrence of the left limit Xpom (t ~). We obtain: 




(Pp Ppf^t , Up PpTp UpC^t ) 

A P {pr,Ur- PrTr) Pp°G‘^ = ’ 


(Pp^ut, Mpout ) = rG(p:j"*)” ,ry^p 

y/y {pr,Ur —PrTr) Otherwise. 

^^g(pOUt)i„ _^^p 

The correctness of these constraints is stated in the following result (see Appen¬ 
dix]^ for a detailed proof): 

Theorem 3. The affine germs of the stationary solutions of the 5-discretization of 


the fluid dynamics are precisely the solutions of System (12) such that Pp,Pq ^ 0 
(pGV, qG Q). 


Since the expressions at the right hand side of the constraints of System (12) 


involve minima of linear terms, these expressions can be interpreted as fractional 


functions over the tropical semifield G. In this way, System (12) can be thought of 


as a set of tropical polynomial constraints (or more precisely, rational constraints). 

The solutions of tropical polynomial systems is a topic of current interest, owing 
to its relations with fundamental algorithmic issues concerning classical polyno¬ 
mial system solving over the reals. Here, we describe a simple method to solve 


System (12), which is akin to policy search in stochastic control. Observe that Sys¬ 


tem (12) corresponds to a fixpoint equation (p, u) = f{p,u), where the function / 


can be expressed as the infimum /\^ f^ of finitely many linear (affine) maps f^. In 
more details, every function f^ is obtained by selecting one term for each minimum 


operation /\ occurring in the constraints (for instance, in (12c), we select one term 
(pp,Up — PpTp) with p G q™). For every selection tt, we can solve the associated 
linear system (p, u) = f^{p, u), and under some structural assumptions on the Petri 
net, the solution {p^,u'^) is unique. If f'^{p'^,u'^) = f{p^,u"), i.e. in every con¬ 
straint, the term we selected is smaller than or equal to the other terms appearing 
in the minimum, then [p^,u'") forms a solution of System ( |12[ ) associated with the 
selection tt. Otherwise, the selection tt does not lead to any solution. Iterating 


this technique over the set of selections provides all the solutions of System (12). 


Every iteration can be done in polynomial time. However, since there is an expo¬ 
nential number of possible selections, the overall time complexity of the method is 
exponential in the size of the Petri net. 
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(13b) {p^^u^) = 


5 . Application to the emergency call center 

We now apply the results of Section]^ to determine the stationary regimes of the 
fluid dynamics associated with our timed Petri net model of emergency call center. 
As in Section [3(4l we consider the subsystem reduced to the variables zi, and zq. 
The corresponding system of constraints over the germ variables {ui,pi), (u^jp^) 
and {uq,pq) is given by: 

(13a) (pi, Ml) = (p5 + TTurPl + TTadvPl, 

Nl + (m 5 - P5Ttr) + 7rur(Mi - PiT^^) + 7radv(Ml - PlTadv)) 

' {P5,N2 + Ms - P5(Ttr + Text)) dv 7rext(Pl, Ml - plText) if P6 = 0 

^T'ext(pi;Mi Pi Text) if P6 ^ 0 

(13c) (P6,M6) = {p 6 ,N 2 - P5(Ttr +T'xt) + (m6 “ Pst'J) A 7rur(pi,Mi - PiTur) 

To solve this system, it is convenient to introduce the following quantity 
T := Text (Text “t” Ttr) T TarTur T TadvTadv 5 

which represents the average time of treatment of a call at level 1 of the model. 
Note that we exclude the trivial case where pi = 0 (and subsequently ps = pe = 0), 
since it cannot occur unless the quantity Ni is null. 

The p-part of (13a) and (13c) show that 

P5 ~ TgxtPl ) 0 ^ Ps ^ '^urPl ■ 

We start by considering the case where pg = 0. Since pi > 0, the minimum 
in (13c I is necessarily attained by the left term. From this, we deduce 

N2 


Pi = 


Text (Ttr + t'xi) ' 


As (psjMs) ^ Text(pi,Mi - piText) (by ( |13b[ )) and ps = TextPi, the inequality Ms < 
Text(Mi — piText) holds. Using the M-part of ( |13a ), we can show that this amounts 
to the inequality 


N2 ^ Text (Ttr H 

— < n := - 


Text) 


We now assume that pg > 0. The fact that Mg = Text(Mi — PiText) (by (13b)) 
leads to the identity 

Ni 

Pi = — ■ 


It remains to distinguish the subcases corresponding to the minimum in (13c). 
• Suppose that the minimum is attained by the left term. We deduce that: 

N2 Ni Text (Ttr + T^xt) ^2 - Nin 
P 6 — — -=-- • 

T T T T 

ur ' ur ' ur 

Since 0 < pg ^ TurPi, we also derive: 

No 






• If the minimum is reached by the right term, then we have pg = TnrPi, or 
equivalently pg = Tur=^. Moreover, we necessarily have uq ^ N 2 — ps(Ttr + T'xt) + 
(Mg — pqtI^j.), which provides ^ ^ r 2 - Note that the latter inequality is strict as 
soon as the minimum in (I3c I is attained by the right term only. 
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Table 1. The normalized throughputs pi , ps and pe as piecewise 
linear functions of N2 / Ni . 


0 < N2/N1 < ri ri < N2/N1 < ra ra ^ N2/N1 


Pi/P* 

Pb/P* 


f N 2 

7rext(Ttr + r'^t) ^1 

f iVa 

nr + r^xt ^1 


P 6 /p* 0 


1 


'^ext 


r 



( 



1 


'^ext 


TTur 


To summarize, we report the possible values of the throughputs pi, ps and pg 
in Table in the stationary regimes. We normalize these values by a quantity p* 
which corresponds to the throughput (of transition qi) in an “ideal” call center 
which involves as many level 2 operators as necessary, i.e. iVa = +oo. Then, the 
throughput p* is given by iVi/f, where f := 7rext{Text + nr) + ’^ur'^ur ^adv^adv 
represents the average time of treatment at level 1. 

As shown in Table the ratios pi /p*, ps/ p* and pg /p* are piecewise linear 
functions of the ratio N 2 /N 1 . The non-differentiability points are given by: 

^ _ 7^'ext(nr + nxt) ^ '^ext(nr + nxt) A '^urT'ur 

f f 

They separate three phases: 

(i) when N 2 /N 1 is strictly smaller than ri, the number of level 2 operators is 
so small that some extremely urgent calls cannot be handled, and no urgent call is 
handled. This is why the throughput of the latter calls at level 2 is null. Also, level 1 
operators are slowed down by the congestion of level 2, since, in the treatment of an 
extremely urgent call, a level 1 operator cannot be released until the call is handled 
by a level 2 operator. 

(ii) when N 2 /N 1 is between ri and r 2 , there are enough level 2 operators to 
handle all the extremely urgent calls, which is why the throughput ps is equal to 
Pi multiplied by the proportion TText of extremely urgent calls. As a consequence, 
level 2 is no longer slowing down level 1 (the throughput pi reaches its maximal 
value p*). However, the throughput of urgent calls at level 2 is still limited because 
N 2 is not sufficiently large. 

(iii) if N 2 /N 1 is larger than r 2 , the three throughputs reach their maximal values. 
This means that level 2 is sufficiently well-staffed w.r.t. level 1. 

This analysis provides a qualitative method to determine an optimal dimension¬ 
ing of the system in stationary regimes. Given a fixed A^i, the number N 2 of level 
2 operators should be taken to be the minimal integer such that N 2 /N 1 ^ r 2 . This 
ensures that the level 2 properly handles the calls transmitted by the level 1 (all 
calls are treated). Then, A^i should be the minimal integer such that pi = ^ 
dominates the arrival rate of calls. 
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Figure 3. Comparison of the throughputs of the non-fluid simu¬ 
lations with the theoretical throughputs (fluid model). The three 
phases are identified by two vertical lines. 


6. Experiments 

We finally compare the analytical results of Section obtained in the fluid 
setting, with the asymptotic throughputs of the Petri net provided by simulations. 


6.1. Asymptotic behavior of the non-fluid dynamics. We have implemented 

the (5-discretization of the non-fluid dynamics (Equations (§-( 111 ), 0 and §) since 
this setting is the closest to reality. Recall that, in this case, tokens are routed 
towards transitions ( 72 , 93 and (74 randomly according to a constant probability 
distribution. We assume that holding times are given by integer numbers of seconds, 
so that we take (5 = Is. In this way, we compute the quantities z^{t) and 

z^{t) by induction on t S N using the equations describing the dynamics. In the 
simulations, we choose holding times and probabilities which are representative of 
the urgency of calls. 

Figure compares the limits when t —>■ -l-cc of the throughputs zi{t)/t, z^{t)/t, 
ze{t)/t of the “real” system, with the throughputs pi, ps and pe of the stationary 
solutions which have been determined in Section The latter are simply computed 
using the analytical formulae of Table [T] We estimate the limits of the throughput 
Zi{t)jt by evaluating the latter quantity for t = 10® s. As shown in Figure]^ these 
estimations confirm the existence of three phases, as described in the previous 
section. The convergence of Zi{t)/t towards the throughputs pi is mostly reached 
in the two extreme phases. In the intermediate phase, the difference between the 
limit of Zi{t)/t and the throughput pi is more important. This originates from the 
stochastic nature of the routing, which causes more variations in the realization of 
the minima in the Zi(t): the throughput of qe increases and the throughputs of qi 
and (75 decrease. 

6.2. Asymptotic behavior of the fluid dynamics. We have also simulated the 
discrete-time fluid dynamics (using Equations §^§ and §-§)■ All simulations 
have been computed with exact rationale in Q. 
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(c) (d) 

Figure 4. Comparison of the fluid dymamics with the stationary 
regimes. Error ratios \zi{t)/t — Pi\/Pi are plotted in log-log scale, 
respectively in blue, red and green when i = 1, 5, 6. 


In most cases, we observe that the corresponding asymptotic throughputs con¬ 
verge to the throughputs of the stationary solutions. This is illustrated in Fig¬ 
ures |3a),|4];b) and|^c), which are obtained using the same set of holding times, 
and by varying the ratio N 2 /N 1 (lower, intermediate and upper phase respectively). 

However, there are also cases in which the convergence does not hold. In the 
experiments we have made, this happens only in the lower phase and in the in- 

Id), 

Sb). 

Such cases suggest the existence of other kinds of stationary regimes of the dynam¬ 
ics, in which the system oscillates between different phases. We remark that the 
non-convergence appears to be related to the existence of arithmetical relationships 
between the holding times of places. An interpretation lies in the fact that, if cycle 
times are not coprime in the system, phenomena of synchronization may lead to 
recurrent slow-down of extreme urgent calls by urgent calls in the two lower phases, 
which could lower the throughput of the system. 


termediate phase, that is, when N 2 /N 1 < r 2 . This is illustrated in Figure 
in which we have increased by one unit of time in comparison to Figure 
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7. Concluding remarks 

We have shown that timed Petri nets with free choice and priority routing can be 
analyzed by means of tropical geometry. This allows us to identify the congestion 
phases in the fluid version of the dynamics of the Petri net. We have applied this 
method to a model of emergency call center. Numerical experiments indicate that 
these theoretical results are representative of the real dynamics. 

In future work, we aim at comparing the behaviors of the fluid deterministic 
model and of the discrete stochastic one. We also plan to study uniqueness condi¬ 
tions of the stationary regimes, and conditions under which convergence of the fluid 
dynamics to a stationary regime can be shown. We will refine our Petri net model 
of emergency call center to take care of the heterogeneous nature of level 2 (calls to 
police and firemen require different instruction times). To this end, it will be helpful 
to implement an analysis tool determining automatically the stationary regimes of 
a timed Petri net given in input. Finally, we plan to analyze the treatment times 
of the system, on top of the throughputs. 
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Appendix A. Proof of Theorem [T] 


Lemma 1. Suppose that all the holding times Tp (p G V) are positive, and consider 
the subpart of the execution trace formed by the transitions fired at the instant t, 
i. e.: 



(with d > 0 unless t = 0, and d' > 0). Then the following two properties hold: 


P+, 


(i) for allp G Vpriority, no transition —> can occur before a transition in (|14|); 
(a) any pair of consecutive of consecutive transitions ^on be switched 

in (14) without changing the states occurring after, provided that [qi,qi+i) 
is not equal to {pff‘*,pf!^*) for some p G Vpriority 


Proof. (i) Suppose that a transition occurs before in (14), i.e. we have 


P- 


pt: 


a subsequence of the form a —> a'... —>. As all the holding times are positive, 
all the tokens consumed by p™* are already present in the state a. In other words, 

. . pT' ■ ■ ■ 

the transition a —> ... is valid in the semantics. This contradicts the priority rule, 
(ii) Consider a pair of consecutive transitions 


such that {qi,qi+i) As discussed in the previous case, since the 

holding times are positive, the transition qi+i does not consume tokens produced by 
the transition qi. Besides, there is no priority rule between qi and qi+i. Therefore, 
qi+i can be fired before qi, and the sequence cr* > leads to the same state 

(7®+^. It follows that all the subsequent states remain identical. □ 


Let us now prove Theorem It is useful to extract the part of the execution 
trace leading to the state cr* of the Petri net at the instant t. It has one of the 
following two forms: 



on 


(16) 


d qi 92 


d‘ 


a 
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depending on whether some transitions q € Q are fired at the instant t or not. In 
both cases, d', d" > 0 and d > 0 unless t = 0, and the durations of the time-elapsing 
transitions occurring before the state cr* in the trace sum up to t. 

In this context, Zq{t) counts the number of transitions occurring before cr* in the 
trace, while Xp{t) is given by the sum of Mp and the number of transitions q S p™ 
occurring before a*. The constraint @ is therefore trivially satisfied by definition 
of Xp(t) and the Zq{t). 

Consider p G T^confiict- Since we use the earliest behavior semantics and p is free 
choice, any token of the initial marking is consumed at the instant 0, and any token 
brought by an upstream transition q' G at the instant s ^ 0 is consumed at the 
instant s + Tp. As a consequence, we can build a bijection which maps each initial 
token with the transition —^ which consumes it at the instant 0, and any transition 

occurring at the instant s — Tp with the transition —^ which consumes at the 
instant s the token brought by q' to place p. We deduce that 

Mp+ Zq>{t-Tp)= Zq{t). 

qf^pin 

Using the constraint this yields to Xp{t — Tp) = J2qep°'‘^ 

Now, let us take q G Qsync- Consider p G Recall that, by definition of Qsync, 
q is the only downstream transition of place p if q. Therefore, every transition —^ 
arising at the instant s consumes a token from place p. This token is either a initial 
token from Mp, or a token brought by a transition q' G p™ fired before the instant 
s — Tp (included). Therefore, we have Zg(t) ^ Xp(t — Tp). In fact, Xp(t — Tp) — Zq(t) is 
equal to the number of tokens with age greater than or equal to Tp located in place 
p in the state cr*. At the instant t + e with 0 < e < d', the age of these tokens will 
be strictly greater than Tp. Therefore, if Xp{t — Tp) — Zq{t) > 0 for all p G g*", the 
transition q can be fired at the instant t + e. But this is impossible in the earliest 
behavior semantics, since the places p G q'’^ are not allowed to contain tokens with 
age strictly greater than Tp while their downstream transition q can be fired. We 
deduce that Xp(t — Tp) = Zqit) for some p G g“. This proves ([^. 

Finally, consider p £ 'Ppriority Using similar arguments as the ones used in the 
previous case, we can show that Zpo^t (t) ^ Xr{t — Tr) for all r G (p™*)™, r ^ p, and 
Zpoat{t) ^ Xr(t — Tj.) for all r G (plU*)™, r ^ p. Besides, we have {t)+ZpMt)^ 
Xp(t — Tp), since every firing of the transition or p™* at the instant s consumes 
a token of Mp or a token brought by an upstream transition of p before the instant 
s — Tp. In consequence, as the function Zpout is non-decreasing, we obtain: 


Zpout (t) + Zpo^t {t ) ^ Xp{t — Tp) . 

In order to prove that § is satisfied, we distinguish two cases depending on the 
form of the trace: 

(i) if the trace is of the form (15), then, by Lemma[^ we can rewrite the subpart 
of the trace as follows: 




times k_ times 


where fc+,fc_ ^ 0, and where p™* and p°“* do not appear in the g'. Then, the 
quantity Xp{t — Tp) Zp<^t ZpOMt {t ) corresponds to the number of tokens with 
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age greater than or equal to Tp in the intermediary state a. If it is positive, and if 
Xr{t — Tr) — Zpo^t (t) > 0 for all r G (p™*)'" such that r ^ p, then the transition p™* 
can be fired right after the state a. This contradicts the priority rule if k- > 0. If 
k_ = 0 , we can fire at the instant t + e {0 < £ < d'), which contradicts the 
definition of the earliest behavior semantics (all the upstream place of contains 
a token older than allowed). 

(ii) if the trace is of the form (16), then Zpout(t“) = Zpout(t). In this case, the 
quantity Xp(t — Tp) — Zp.^t(t) — Zpaat{t~) represents the number of tokens with age 

■rt 


IpJ yv J 

greater than or equal to Tp at place p in the state ct*. If Xp{t — Tp) — Zp<^t[t) — 
Zpout(t“) > 0 and Xr{t — Tr) — Zp<^t{t) > 0 for all r G such that r p, then 

the transition p™* can be fired at the instant t + e with 0 < e < d'. This is again 
a contradiction with the earliest behavior semantics. 

In both cases, we have Xp{t—Tp)—Zp^t (f)—z^out (t~) = 0 or Xr(t—Tr)—Zp^t (t) = 0 




for some r G such that r ^ p. We deduce that the constraint holds. 

Now assume that Xp{t — Tp) — Zp^t(t) — z^out (t) > 0 and Xr{t — Tr) — Zpout (t) > 0 
for all r G such that r ^ p. These quantities correspond to the number of 

tokens in places p and r with age greater than or equal to Tp and Tr respectively, 
in the state cr*. Thus, the transition p°“* is activated at the instant t + e for 
all e > 0 sufficiently small. Note that Xp{t — Tp) — Zpmit(t) — Zpout(t“) > 0 as 
ZpOUt {t ) < Zp out (t). Thus, there exists a place r' G (p™*)'" with r ^ p, such that 
Xr'(t — Tr') = ZpOUt (t). In Other words, place r' does not contain any token with 
age greater than or equal to r^'. Given e > 0 sufficiently small, this is still true at 
the instant < + £, so that the transition p™*^ cannot be fired at the instant t + e. 
Therefore, we are allowed to fire the transition p°“* at the instant t + e, which is 
a contradiction with the definition of the earliest behavior semantics. As a result, 
0 is satisfied. □ 


Appendix B. Proof of Theorem [3] 


We denote the lexicographic order over by and we use the notation x < y 
when X ^ y and x ^ y. 

We first remark that given two functions / and g of affine germs (p, u) and 
{p', u'), the function t i-A f(t)Ag{t) belongs to the affine germ given by the minimum 
{p,u) A {p',u') taken in the lexicographic order. Besides, if r G (5N, the function 
1 1 —>■ f{t — t) belongs to the affine germ (p, u — Tp). Finally, for all A G M, the affine 
germ of the map t i-A Xf{t) is equal to A(p, u) = (Ap, Xu). 

Now, suppose that Xp and Zq are stationary solutions of the d-discretization of 
the fluid dynamics, and let (pp, Up) and (p^, Uq) be the respective germs, for p G P 


and q G Q. Since the functions Xp anui .oq 


and Zq satisfy the constraints ([^, (|^, (|7|, ( |^ for 
all t G 5N, we deduce from the previous properties that the constraints (12a|-(12d) 
are satisfied. Besides, given p G T^priority, 


ensures that for all t G SN, we have: 


Zp^ut(t) = [xp{t - Tp) - Zpaut{t - S)) A 


A 


Xr{t — Tr) . 


' , rytp 
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Consequently, we obtain 
(17) 

) (Pp Ppojit , Up PpTp Upojit “t" Ppo^t5^ A l\^ ^Pr-i Pr'^r) ■ 

'’ 6 (p+“*)‘" . ■'■Ap 


If Ppout = 0, this amounts to the constraint given in (12 51. Now consider the case 


where Ppout > 0. Using (12d), we know that 

(Pp^ut , MpOUt ) (Pp Pp^ut , Up PpTp ZlpOUt ) . 


and thus 


(PpT^* 5 ^p5^^) ^ (Pp Ppyi^t 5 ^p Pp'^p Up^ut) . 


Since d > 0, it follows that 

(Pp^ut , Upout) ^ (pp — Ppout , Up — PpTp — Upout + PpC^tS) . 

We conclude that the constraint 0 is equivalent to ( |12^[ ) when ppout > 0. 

Conversely, let (pp,Mp) and {pq,Uq) be solutions of System ( [T^ . We define Xp 
and Zq as the functions given by Xp(t) = Up + Ppkd and Zq{t) = Uq + pqkd for all 
t G [k5,{k + l)d) and fc S N. The constraints (12a)"(12d) ensure that (§, (§, 
(0, @ hold for all t G (5N. Since all the holding times Tp belong to dN and the 
functions Xp and Zq are constant on the intervals of the form [kS, (fc+l)(5), we deduce 
that these constraints 0, 0, 0, 0 actually hold for all t in such intervals, and 
so for all t ^ 0. Moreover, as previously shown, the constraint given in (12 5) 
is equivalent to ( [l7| ), since ( |12d ) is satisfied. This proves that the constraint ([8) 
holds for all t G dN. It remains to show that the latter constraint is satisfied 
when t G (k6, {k + l)d). First observe that Zpo^tikS) ^ ArG(p°"*)‘" r^p^rikS — Tr) 
ensures that Zp.^t{t) ^ ArG(p5"‘)“ r^p^r{t — t^)- Besides, by 0, we know that 
Zpout(t) ^ Xp{t — Tp) — Zpaat(t). We now distinguish two cases: 

(i) if we have Zp<^t[k6) = ArG(p°“‘)‘" r^p^^ikS — Tr), then straightforwardly, 
Zpout (t) = ArG(p!j"*)'“ , r^p ~ U)- 

(ii) if Zpout(fcd) = Xp(k5 — Tp) — Zpont((k — l)i5), we obtain: 

Zpa^t(k5) ^ Xp(k5 — Tp) — ZpaMt(kS) ^ Zpa^t(k5 ), 

where the first inequality comes from the fact that Zpo^t is non-decreasing, and the 
second inequality from 0. We deduce that Zpout(fc(5) = Xp{kS — Tp) — Zpant(^k6). 
Hence, we get Zpout(f) = Xp{t — Tp) — Zpout(t). 

As a consequence, in both cases, we have proved that Zpout {t) is the minimum 
between Zpout(t) — Zpout(t) and ArG(p°"‘)‘“ r^p^rit — Tr). This shows that 0 holds 
for all t > 0. □ 
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